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Radiative cooling is central to a wide range of astrophysical problems. Despite its 
importance, cooling rates are generally computed using very restrictive assumptions, 
such as collisional ionization equilibrium and solar relative abundances. We simul- 
taneously relax both assumptions and investigate the effects of photo-ionization of 
heavy elements by the meta-galactic UV/X-ray background and of variations in rela- 
tive abundances on the cooling rates of optically thin gas in ionization equilibrium. We 
find that photo-ionization by the meta-galactic background radiation reduces the net 
cooling rates by up to an order of magnitude for gas densities and temperatures typ- 
ical of the shock-heated intcrgalactic medium and proto-galaxies (10^K<r< 10^ K, 
pI (p) ^ 100). In addition, photo-ionization changes the relative contributions of dif- 
ferent elements to the cooling rates. We conclude that photo-ionization by the ionizing 
background and heavy elements both need to be taken into account in order for the 
cooling rates to be correct to order of magnitude. Moreover, if the rates need to be 
known to better than a factor of a few, then departures of the relative abundances 
from solar need to be taken into account. We propose a method to compute cool- 
ing rates on an element-by-element basis by interpolating pre-computed tables that 
take photo-ionization into account. We provide such tables for a popular model of 
the evolving UV/X-ray background radiation, computed using the photo-ionization 
package cloudy. 
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1 INTRODUCTION 

Dissipation of energy via radiative cooling plays a central 
role in many different astrophysical contexts. In general 
the cooling rate depends on a large number of parameters, 
such as the gas density, temperature, chemical composi- 
tion, ionization balance, and the radiation field. In the 
absence of radiation, however, the equilibrium ionization 
balance depends only on the temperature. In that case 
the coohng rate in the low density regime, which is dom- 
inated by collisional processes, is simply proportional to 
the gas density squared, for a given composition. Thus, 
the coohng rates for a plasma in collisional ionization 
equilibrium (CIE) can be conveniently tabulated as a 
function of the tem perature and composition (metallic- 
ity) of the gas (e.g., ICox fc Tucker! [l969l : [Raymond et al.l 
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!2008! ) 



12001!; 

such 



and 
tables are widely used for a large variety of problems. 

Although it is convenient to ignore radiation when cal- 
culating cooling rates, radiation is generally important for 
the thermal and ionization state of astrophysical plasmas. 
For example, [Efstathiou (1992,) investigated the effect of the 
extragalactic UV background on the cooling rates for gas 
of primordial composition (in practice this means a pure 
H/He plasma) and found that including photo-ionization 
can suppress the cooling rates of gas in the temperature 
range T ~ 10* — 10^ K by a large factor. Although the ef- 
fects of radiation are often taken into account for gas of 
primordial composition, photo-ionization of heavy elements 
is usually ignored in the calculatio n of cooling rates (but see 
!Leonard|[l99i : !Benson et aLl!2002! ). 
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In this paper we will investigate the dependence of cool- 
ing rates of gas enriched with metals on the presence of 
ionizing radiation, focusing on the temperature range T '^ 
10^ — 10* K and optically thin plasmas. We will show that, 
as is the case for gas of primordial composition (|Efstathioij 
1 1992 ). photo- ionization can suppress the metallic cooling 
rates by a large factor. Moreover, the suppression of the 
cooling rate is significant up to much higher temperatures 
than for the primordial case. 

We will also investigate the relative contributions of var- 
ious elements to the cooling rates. If the relative abundances 
are similar to solar, then oxygen, neon, and iron dominate 
the cooling in the temperature range T ~ 10* — lO'^ K. How- 
ever, we will show that the relative contributions of difi'erent 
elements to the cooling rate are sensitive to the presence of 
ionizing radiation. 

Although we will illustrate the results using densities 
and radiation fields that are relevant for studies of galaxy 
formation and the intergalactic medium (IGM), the conclu- 
sion that photo-ionization significantly reduces the cooling 
rates of enriched gas is valid for a large variety of astrophys- 
ical problems. For example, for T ~ lO'^ K and T ~ 10*^ K 
the reduction of metal-line cooling rates is significant as 
long as the dimensionless ionization parameteijj U > 10~^ 
and [/>10~^, respectively. We will focus on the temper- 
ature range 10** — 10* K because gas in this temperature 
range is usually optically thin and because the effects of 
photo-ionization are generally unimportant at higher tem- 
peratures. 

Tables containing cooling rates and several other use- 
ful quantities as a function of density, temperature, red- 
shift, and composition, appropriate for gas exposed to 
the model s for the evolving meta- galactic UV/X-ray back- 
ground of iHaardt fc Madaul (|200in are available on the fol- 
lowing web site: http://www.strw.leidenuniv.nl/WSS08/. 
The web site also contains a number of videos that illustrate 
the dependence of the cooling rates on various parameters. 

This paper is organized as follows. In Section [5] we 
present our method for calculating element-by-element cool- 
ing rates including photo-ionization and we compare the 
limiting case of CIE to results taken from the literature. 
Section[3]shows how metals and ionizing radiation affect the 
cooling rates. Section |4] demonstrates the importance for the 
low-redshift shock-heated IGM, which is thought to contain 
most of the baryons. In this section we also illustrate the ef- 
fect of changing the intensity and spectral shape of the ion- 
izing radiation. We investigate the effect of photo-ionization 
on the relative contributions of individual elements in Sec- 
tion [5] and we summarize and discuss our conclusions in 
Section |B1 

Throughou t this paper w e use the cosmological pa- 
rameters from iKomatsu et al.l l|2008l ): {Q,in,i^A,^h,h) = 
(0.279,0.721,0.0462,0.701) and a primordial helium mass 
fraction Xhc ~ 0.248. Densities will be expressed both as 
proper hydrogen number densities nn and density contrasts 
5 = pb/ (Pb) — 1, where (pb) is the cosmic mean baryon 
density. The two are related by 



Table 1. Default CLOUDY solar abundances 



nn 



2 METHOD 



Element 


Ui/nii 


Mass Fraction 


H 


1 


0.7065 


He 


0.1 


0.2806 


C 


2.46 X 10-" 


2.07 X 10-=" 


N 


8.51 X 10-5 


8.36 X 10-* 


O 


4.90 X 10-'' 


5.49 X 10-^ 


Ne 


1.00 X 10-4 


1.41 X 10-^ 


Mg 


3.47 X 10-5 


5.91 X 10-* 


Si 


3.47 X 10-5 


6.83 X 10-* 


S 


1.86 X 10-5 


4.09 X 10-* 


Ca 


2.29 X lO-is 


6.44 X 10-5 


Fe 


2.82 X 10-5 


1.10 X 10-^ 


X 10"'^ cm 


-^ ri +S)C\-i 


-.i«r^« V 



(1) 



All radiative cooling and heating rates were computed by 
running large grids of photo-ionization models using the 
publicly available photo-ionization pa ckage CLOUD'vn (ver - 
sion 07.02 of the code last described bv lFerland et al.lll998l ). 
CLOUDY contains most of the atomic processes that are 
thought to be important in the temperature range of inter- 
est here (T ^ 10* — 10* K) and the reader is referred to the 
online documentation for details about the atomic physics 
and data used. 

The gas was exposed to the cosmic microwave back- 
ground radiation (CMB) and the iHaardt fc Madaul (|200ll . 
hereafter HMOl) modeo for the UV/X-ray background ra- 
diation from galaxies (assuming a 10 percent escape fraction 
for H-ionizing photons) and quasars. We assumed the gas to 
be dust-free, optically thin and in ionization equilibrium. 
We discuss the limitations and the effects of the last two 
assumptions in section [S] All cooling rates are tabulated as 
a function of lognn (total hydrogen number density), logT 
(temperature), and z (redshift). 

When studying the effect of changes in the relative 
abundances of elements, we compute the cooling rates 
on an element-by-element basis. The cooling rate Ai (in 
ergs-'^ cm-"^) due to element i, where element i is heav- 
ier than helium, is defined as the difference between the 
cooling rate computed using all elements (assuming solar 
abundances) and the cooling rate computed after setting 
the abundance of element i to zero, while keeping all other 
abundances (i.e., number densities relative to H) the same. 
This is a valid approximation provided element i does not 
contribute significantly to the free electron density, which 
is the case for all elements heavier than helium and if the 
metallicity Z < Zq . 

The combined contributions from hydrogen and helium 
are computed by interpolating in the four dimensions log riH , 
logT, z, and nuc/nu from tables of CLOUDY models that 
contain only H and He. 

Thus, the total net cooling rate. 



^ U = 3>H/(riHc), where <I>h is the flux of hydrogen ionizing pho- 
tons (i.e., photons per unit area and time), nn is the total hydro- 
gen number density and c is the speed of light. 



^ http://www.nublado.org/ 

•^ http : //pitto . mib . inf n . it/~haardt/ref model . html 
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A = Ah,ho+ ^ Ai, 

i>He 

could be obtained from 



A = At 



V- 



(ni/nn)© 



A»,e 



(2) 



(3) 



where (ni/nn)© is the solar abundance of element i (the 
default CLOUDY solar abundances are given in Table [TJ and 
Ai,0 is the contribution of heavy element i to the radiative 
cooling rate for solar abundances, which we have tabulated 
as a function of lognn, logT, and z. Note that we use A to 
denote the cooling rate per unit volume (ergs~^ cm^'^). 

We can, however, do better than equation (|3]) by tak- 
ing the dependence of the free electron density on He/H, 
into account (in equation @, the electron density is im- 
plicitly assumed to be that corresponding to solar abun- 
dances - nHe/nn = 0.1). Since cooling rates due to metals 
are dominated by collisions between ions and free electrons, 
Ai scales as the product of the free electron and ion densities, 
Ai oc Ticrii. Hence, 



A = Ah, He + ^ A, 

i>Hc 

= Ah, He + > Ai,^y 

^-^ [ne/nnJQ 



n-c,© 
Uc/nu 






10'* 



/H] 



(4) 



where lo'*'^' = (ni/nH)/(ni/nH)© and we used the fact that 
riH = nn,© (since we tabulate as a function of nn). While 
(rie/nn)© is obtained by interpolating the solar abundance 
table for n-o/riH in {n-a,T,z), (ne/nn) must be obtained by 
interpolating in ((r^Hc/nH),?^H,^, z). Note that we tabulate 
the electron density in the absence of metals. This is a valid 
approximation given that heavy elements only contribute 
significantly to the free electron density for Z ^ Zq. 

In practice we tabulate logAi,©/nH over the range 
log[rzH (cm-^)] = -8.0, -7.9, . . ., 0.0 and log[r (K)] = 2.00, 
2.02, . . ., 9.00. In addition, the quantities logAH,Hc/nH) 
rLc/n-a and the mean particle mass are all tabulated as a 
function of density and temperature for each of the values 
riHe/riH = 0.0787, 0.0830, 0.0876, 0.0922, 0.970, 0.102, 0.107 
which correspond to mass fractions Xhc/(A'h + Xhc) = 
0.238, 0.248, 0.258, 0.268, 0.278, 0.288, 0.298. Finally, for 
each value of the ratio tiho/wh we tabulate the temperature 
as a function of density and internal energy per unit mass 
to enable simulation codes that parametrize thermal energy 
in terms of the latter quantity to use the cooling tables. 
We have computed these tables for each of the 49 redshifts 
spanning z = — 9 for which the HMOl models are defined. 
All tables are in HDF5 format and together they result in 
upwards of 232 MB of storage. Reducing the resolution in 
riH and T by a factor of two does reduce the accuracy of 
interpolated rates significantljQ, but reduces the storage re- 
quirements to 61 MB. 

We find that for metallicities Z < Zq equation (|4]) 



* Reducing the temperature resolution by a factor of two roughly 
doubles the interpolation errors near the thermal equilibrium so- 
lution, but in this regime the cooling times are in any case effec- 
tively infinite. 




Figure 1. Comparison of normalized CIE cooling rates 
for various studie s . All curves use the solar abundances of 
iGnat fc Sternberg 112003), which differ somewhat from our de- 
fault solar abundances. Sho wn are CLOUDY version 07.02 (solid), 
[Sutherland fc Dopital lll993l) fd ashed). MAPPINGS ni (dot-dashed), 
and iGnat fc Sternberj 12003) (dot-dot-dot-dashed). MAPPINGS 
gives significantly higher cooling rates for T ~ 10^ K, but the 
differences are typically smaller than a factor of two. Note that 
for this comparison the cooling rates were divided by nnfte, but 
that we divide by n^ in our tables and in figure [6] The upturn in 
the normalized cooling rates below 10** K is caused by the sharp 
decrease in rio with decreasing temperature. 



closely matches the true cooling rate (i.e., including all ele- 
ments) with very good accuracj|f| if the following 11 elements 
are included: H, He, C, N, O, Ne, Mg, Si, S, Ca, and Fe. 

For redshift z = and over the full range of densi- 
ties and temperatures, the median relative errors in the 
absolute net cooling rates are 0.33%, 1.6%, and 6.1% for 
Z — O.IZq, Zq, and WZq, respectively (we have scaled 
nHc/riH with metallicity) . For higher redshifts the median 
errors are smaller than for z — because Compton cooling 
off the CMB, which is modelled accurately, becomes increas- 
ingly important. Hence, even for metallicities as extreme as 
10 times solar, using equation Q and including 11 elements 
gives errors of only a few percent. As we shall see below, 
this is much smaller than the differences between different 
photo-ionization codes. Using equation (|3]) rather than Q 
gives similar errors for low metallicities, but the median er- 
rors are a factor 2-3 higher for Z > lO"'^ Zq. 

Excluding temperatures within 0.1 dex of the thermal 
equilibrium solution (where the relative errors in the net 



^ While the error in the total cooling rates is small, the relative 
errors in the contribution of individual elements can be larger if 
the element contributes negligibly to the total cooling rate (see 
e.g., the noise in the Fe contour at T < 10^ K in the left-hand 
panel of figure (BJl. This is because the tables for individual ele- 
ments were computed by taking the difference between the total 
cooling rate for solar abundances including and excluding the el- 
ement. Hence, if the element does not contribute significantly to 
the cooling rate for a particular density, temperature, and red- 
shift, then its contribution will be computed as the difference 
between two nearly equal numbers. Because this can only occur 
if the contribution is negligible, we do not consider this to be a 
problem. 
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cooling rate, which is computed as the absolute difference 
between heating and cooling, become large because we are 
subtracting two nearly identical numbers), the maximum 
errors in the net cooling rates are 32%, 29%, and 39% for 
Z — O.IZq, Zq, and IOZq, respectively. Thus, even for ex- 
treme metallicities the maximum difference between the es- 
timated and true rates is well within a factor of two. These 
maxima are reached near the 0.1 dex exclusion zone, so the 
maximum errors for temperatures that differ substantially 
from the equilibrium values are much smaller. 

We have also tabulated the combined cooling rates of all 
elements heavier than helium, assuming solar relative abun- 
dances. While scaling all heavy elements simultaneously re- 
duces the accuracy (relative to scaling the contribution of 
each element individually), these tables may be convenient 
in the absence of a complete set of elemental abundances. 
In this case, the total cooling rate becomes: 



A = A 



H.Ho 



+ A 



Z,Q 



Uc/riH 



(nc/nn)© Z@ 



(5) 



Although cooling rates including photo-ionization have 
not yet been tabulated, CIE cooling rates have been pub- 
lished by various authors. In Figure [l] we compare our 
CIE results (i.e., clou dy version 07.02) with those of 
Gnat fc Sternberg! I^2007^ (w ho used cloudy version 06.02), 
Sutherland fc Dopital (|l993) (w ho used mapping s) and also 
with MAPPINGS III (version r: Groves et aLllJOOSl ). The cal- 
culations generally differ in terms of the code, the atomic 
rates, and the solar abundances that were used. In order to 
focus the comparison on c odes and rates, we used the so- 
lar abundances assumed bv lGnat fc Sternberg! (|2007|) when 
computing the cooling rates for all curves in this figurqj. 
Note that these abundances differ somewhat from the solar 
abundances that we use in the rest of this paper. 

The differences in the cooling rates shown in Figure [T] 
are typically smaller than 0.3 dex suggesting that the cool- 
ing rates are known to better than a factor of two, at least 
for CIE. The differences are largest for T ~ 10^ K, with 
MAPPINGS giving higher cooling rates than cloudy. 



3 PHOTO-IONIZATION, METALS, AND 
COOLING RATES 

The cooling rates in photo-ionization equilibrium are not 
proportional to the density squared, as is the case for CIE. 
As a consequence, the cooling curve changes dramatically 
with density. This is illustrated in Fig. [2] which shows nor- 
malized, net cooling, rates (lA/n^j), as a function of tem- 
perature for solar abundances. Different curves correspond 
to different densities, as indicated in the legend. 

The calculation includes the CMB and the HMOl back- 
ground at redshift three. The z — 3 HMOl background 
corresponds to a hydrogen photo-ionization rate ri2 = 



r/lO" 



— 1.1 and an ionization parameter U — 



1.6 X 10 ^/tih. We stress that in the vicinity of ionizing 
sources the radiation field may be much more intense, which 



^ ISutherland &: Dopital 1119931 ) only give the cooling rates as a 
function of th e abundance of iron whi ch we interpolated to the 
corresponding iGnat &: Sternberei 1120071 ) value. 
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Figure 2. Normalized, absolute, net cooling rates (|A/n^|) as 
a function of temperature for solar abundances. Different curves 
correspond to different gas densities, as indicated in the legend. 
The gas in all models is exposed to the z = 3 CMB and the 2 = 3 
HMOl UV/X-ray background. Note that the cooling rates were 
divided by n^, but that we divided by nufie in Fig.lTI 



would enhance its effect on the cooling rates compared to 
the results presented here. 

As the density decreases the gas becomes more highly 
ionized and the cooling peaks due to coUisional excitation of 
various ions disappear. For very low densities and high tem- 
peratures the normalized cooling rates actually increase as 
they become dominated by Compton cooling off the CMB, 
which scales as A oc ricT. For the two lowest densities shown, 
the equilibrium temperature, for which the net cooling rate 
A ^ 0, exceeds 10'' K. In these cases, curves corresponding 
to net heating appear to the left of the thermal equilibrium 
temperature. 

Figure [3] shows how the presence of heavy elements 
and radiation affects the radiative cooling time. Each panel 
shows contours of constant cooling time in the density- 
temperature plane for two models. The top panels illustrate 
well known results, while the bottom panels demonstrate the 
importance of radiation on the cooling due to metal lines. 

Here and throughout the cooling time refers to the ab- 
solute value of the net radiative cooling time at a fixed hy- 
drogen density. 



f-rool — 



T 



\nkT 



dT/dt |Ah 



Ac 



(6) 



For temperatures below the thermal equilibrium tempera- 
ture, fcooi corresponds to a heating timescale. The calcu- 
lation includes the CMB and, optionally, the HMOl back- 
ground at redshift three. 

Note that adiabatic cooling due to the expansion of the 
Universe is not included, but would dominate over radiative 
cooling for sufficiently low densities. For example, the uni- 
versal Hubble expansion (appropriate for densities around 
the cosnnc mean, i.e., nn ^ 10"^ cm~'^ at z = 3) corre- 
sponds to an adiabatic cooling time tcooi.adiab = 1/(2//) ~ 
1.6 X 10^ yr at z = 3. 

The top-left panel shows that for a metal-freqj gas. 



We use riHc/'^H = 0.083 for primordial abundances. 



Cooling of enriched, photo-ionized gas 5 





Z = Za. r,2 = 1.1 
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Figure 3. Contour plots of the net cooling time in years as a function of temperature and total hydrogen number density (note that 
the mean baryon density corresponds to tih ~ 10~^ cm~^ at 2: = 3). For temperatures below the equilibrium temperature, the contours 
indicate negative net cooling times, that is, heating times (see text). The gas in all models is exposed to the 2 = 3 CMB. In addition, 
models for which the legend indicates ri2 = 1.1 use the z = 3 HMOl UV/X-ray background. The figure confirms that both metals and 
photo-ionization affect the cooling rates significantly, and shows their combined effect. Top-left: The effect of photo-ionization for zero 
metallicity; Top-right: The effect of metallicity for a purely collisionally ionized plasma; Bottom-left: The effect of photo-ionization for 
solar metallicity; Bottom-right: The effect of metallicity for a photo-ionized plasma. 



ionizing radiation drastically reduc es the cooling rates for 
10* K < T< lO^K (|Efstathioulll992l ). This happens because 
the peaks in the cooling rate due to collisional excitation 
of neutral hydrogen and singly ionized helium (followed by 
Lya emission) are removed when the gas is photo-ionized. 

In the presence of ionizing radiation (solid curves) there 
are actually two contours for each net cooling time, cor- 
responding to net heating (left of the thermal equilibrium 
asymptote) and net cooling (right of the equilibrium asymp- 
tote). These two contours nearly mergqj near the (density- 
dependent) thermal equilibrium temperature, slightly above 
10* K, where the net cooling time goes to infinity. Below 
this temperature, the heating time is dictated by the photo- 
heating rate. 

For T S> 10^ K photo-ionization has no effect on a 
metal-free gas because the plasma is already fully ionized 
by collisional processes. In this regime the plasma cools pre- 



* The fact that the contours really do merge and then disappear 
at high densities is due to the limitations of our plotting package 
and the finite resolution of our grid. 



dominantly via the emission of Bremsstrahlung and/or in- 
verse Compton scattering of CMB photons. The latter pro- 
cess is the dominant cooling mechanism for most of the 
baryons at high redshift {z > 7). In the plot inverse Comp- 
ton coohng off the CMB dominates the radiative cooling 
rate at low densities and high temperatures, but the corre- 
sponding cooling time exceeds the Hubble time. 

For diffuse, intergalactic gas (i.e., for density contrasts 
5 <^ 10^, corresponding to riH ^ 10"'^ cm"'' at 2; = 
3, much lower than expected in virialized objects, e.g., 
IColes fc Lucchinll2002r ) of primordial composition, radiation 
increases the cooling time at 10^ K by at least an order of 
magnitude and by much more at lower temperatures. Since 
the cooling times in this regime are comparable to the Hub- 
ble time, radiation will have a large effect on the fraction 
of the baryons that are hot. At densities corresponding to 



collapsed objects (5 > 10 or nu > 10 



: 3), the 



increase in the cooling time is generally smaller, although it 
can still easily be an order of magnitude at temperatures as 
high as IC'-'^K. 

The top- right panel shows that heavy elements strongly 
increase the cooling rate of a collisionally ionized plasma 
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for 10'' K < T<10'^K (e.g.. iBoehringer fc HensIe^llQSgl ). 
Comparing the model with primordial abundances (dashed 
contours) to the one assuming solar metallicity (solid con- 
tours), we see that the cooUng times typically differ by about 
an order of magnitude. The presence of metals allows radia- 
tive cooling through coUisional excitation of a large number 
of ions at a variety of temperatures. For T > 10^ K the dif- 
ference is smaller because there are few lines to excite since 
most elements are collisionally ionized to a very high degree. 
Bremsstrahlung is the dominant cooling mechanism at these 
very high temperatures. 

The bottom-left panel demonstrates that ionizing ra- 
diation also strongly reduces the cooling rates when heavy 
elements dominate the cooling. This happens for the same 
reason as in the primordial case. The radiation field ion- 
izes the plasma to a higher degree than it would be in CIE. 
Hence, the ions that are typically collisionally excited are 
not present, reducing the cooling rates. 

Note that a similar thing happens when a colli- 
sionally ionized plasma is cooling more quickly than 
it can recombine. In that case the ionization balance 
shifts out of equilibrium leaving the gas too highly ion- 
ized for its temperature. As for photo-ionization, the 
associated reduction of the number of bound electrons 
with excitation energies low enough to be collis i onally 
excited, reduces the c o oling rate (^e.g.. [Ka fatos 1973; 
Shapiro fc Moore 19761: Schmutzler fc Tscharnuter 1995 : 



Sutherland fc Dopitalll993[ r 

Comparing the models with (solid contours) and with- 
out ionizing radiation (dashed contours), we see that for den- 
sities characteristic of the diffuse IGM {S ^ 10^ correspond- 
ing to riH <^ 10~^ cm~^ at z; = 3), radiation significantly re- 
duces the cooling rates for T < 10® K and that the reduction 
typically exceeds one order of magnitude for r< 3 x 10^ K. 
For densities nn ^ 10~^ cm~^ the cooling rates are only sup- 
pressed substantially for T < 3 x IC* K. However, we stress 
that what matters here is the ionization parameter. Hence, 
a stronger radiation field will affect the cooling rates up to 
higher densities. Gas with a density characteristic of col- 
lapsed objects will typically be close to sources of ionizing 
radiation and may thus be exposed to a radiation field that 
is more intense than the meta-galactic UV background. 

Finally, the bottom-right panel shows that metals 
strongly increase the cooling rates in the presence of an ion- 
izing radiation field, as was the case for a collisionally ionized 
gas (top-right panel). For solar metallicity the cooling rate 
at r < 10^ K is typically an order of magnitude higher than 
for a primordial composition. 

Note that for very low densities (nn ^ 10~^ cm~^) met- 
als increase the equilibrium temperature (towards the left 
side of the region where the net cooling time is above 10^^° 
years) because their presence boosts the photo-heating (via 
oxygen and iroqj) but does not significantly affect the cool- 
ing because it is dominated by Compton cooling off the 
CMB. 




Figure 4. Contours show where the net cooUng time equals 
the Hubble time in the density contrast - temperature plane at 
2 = 0. From top to bottom, the different colours correspond to 
primordial composition (blue), a metallicity of ten percent so- 
lar (red), and solar abundances (black). Dashed contours are for 
coUisional ionization only, while solid contours arc for gas ex- 
posed to the HMOl model for the z = UV/X-ray background. 
Each solid contour comprises two components, corresponding 
to not heating (low temperature component) and net cooling 
(high temperature component), respectively, which merge near 
the (density-dependent) thermal equilibrium temperature. Both 
photo-ionization and metallicity determine whether gas is able to 
cool and it is therefore crucial to take both into account when 
predicting the fraction of the baryons that reside in the WHIM. 



4 EFFECT ON THE WHIM 

The fraction of gas that has been shock-heated to tem- 
peratures of 10^ K < T < 10^ K is currently of great in- 
terest, mainly because this so-called warm-hot intergalac- 
tic medium (WHIM) is hard to detect, yet may contain a 
large fraction of the baryo ns in the low-redshift Universe 
(e.g.. lCen fc Ostrikerlll999r ). There are two reasons why the 
WHIM becomes more important at lower redshift. First, as 
structure formation progresses, larger structures form, lead- 
ing to stronger gravitational accretion shocks and a greater 
fraction of the baryons are heated to temperatures in the 
WHIM range. Second, as the universe expands, the density 
of the diffuse gas decreases a^ij (1 -1- z)^ and the cooling 
time due to coUisional processes (which dominate for z < 7) 
will thus increase as (1 + z)~^. Hence, the cooling time in- 
creases faster than the Hubble time and more and more of 
the shock-heated gas is unable to cool. 

Since the cooling times are of order the Hubble time for 
much of the WHIM, the precise values of the cooling rate 
is particularly important. Because the cooling of the WHIM 
tends to be dominated by line radiation, beca use its density 
is low (5 ~ 10 — 10^; e.g.. lBertone et al.ll2008l ). and because 
the WHIM gas may well be enriched to values of 10% of 
solar or higher, both metals and photo-ionization by the 
UV background may be important. Figure |4] demonstrates 
that this is indeed the case. 

Figure m shows a contour plot of the cooling time in the 



^ The increased He/H ratio also boosts the photo-heating rate, 
but its contribution is smaller than that of the metals. 



^'^ In reality the density will decrease slightly less fast with time 
if the gas is overdense and collapsing. 
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density contrast-temperature plane for redshift z = 0. Each 
contour corresponds to the same net coohng time, namely 
the Hubble time. Dashed contours are for a purely coUi- 
sionally ionized gas, while the solid contours include photo- 
ionization by the z — UV/X-ray background radiation 
from galaxies and quasars, which corresponds to a hydro- 
gen photo-ionization rate ri2 = 0.08 and an ionization pa- 
rameter [/ = 9 X 10~'^/nH. Metallicity increases from the 
top down from zero (i.e., primordial abundances) (blue con- 
tours), to 10% solar (red) to solar (black), as indicated in 
the figure. Gas above a given contour is able to cool, while 
gas below it will remain hot. 

For a primordial composition (blue contours), turning 
on the ionizing radiation raises the density contrast above 
which the gas can cool from 5 < 10 to 5 ~ 10^ for T < lO'^ K. 
For solar metallicity, the UV radiation becomes important 
for T < 10^ K and raises the critical density required for 
cooling within a Hubble time by about an order of magni- 
tude for T<10"'K. Increasing the metallicity from zero to 
solar decreases the critical density contrast by about an or- 
der of magnitude over the full range of WHIM temperatures. 
Clearly, both photo-ionization and metals are important for 
the thermal evolution of the WHIM. 

Figure[5]shows how the intensity and the spectral shape 
of the ionizing radiation affect the results for solar abun- 
dances. The solid contour indicates the density contrasts and 
temperatures for which the cooling time equals the Hubble 
time for the z — HMOl model for the UV/X-ray radiation 
from quasars and galaxies (model QG). It is identical to the 
black, solid contour in Figure |4] The dashed contour shows 
the fcooi = tn contour after we have multiplied the entire 
HMOl radiation field by a factor of 10 (model 10*QG). In- 
creasing the intensity reduces the cooling rates, shifting the 
contour to higher densities. 

The dot-dashed contour corresponds to the quasar only 
HMOl model (model Q), which we have rescaled such as to 
give the same hydrogen ionization rate as model QG. Since 
the average energy of H-ionizing photons is higher for model 
Q than for QG, the difference between the dot-dashed and 
solid contours reflects the effect of changing the spectral 
hardness. Comparing these two contours, we see that, at a 
fixed H ionization rate, harder spectra tend to inhibit the 
cooling more than softer spectra for temperatures above the 
equilibrium value (the tcooi = tn contour moves to higher 
densities for harder spectra). This happens because coUi- 
sional excitation of heavy ions involves electrons with higher 
ionization energies than hydrogen. The extra high energy 
photons will remove many of those electrons, thus reduc- 
ing the cooling rates further. For temperatures below the 
equilibrium value, however, the net cooling time actually 
decreases (the lower-left of the left contour moves to lower 
densities for harder spectra) because the photo-heating rates 
are larger for harder spectra. 

In summary, for densities and temperatures character- 
istic of the WHIM the nature of the ionizing radiation field 
as well as the metallicity of the gas may have a significant 
impact. Accurate cooling rates therefore require a correct 
treatment of the composition of the gas, the spectral hard- 
ness, and the radiation intensity. Simple models may result 
in a poor estimate of the amount of baryons in the WHIM 
phase. 
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Figure 5. Contours show where the net cooling time equals the 
Hubble time in the density contrast - temperature plane for gas 
of solar metallicity at 2 = 0. The solid contour is for the HMOl 
model for the z = UV/X-ray background from galaxies and 
quasars (model QG), the dashed contour is for a radiation field 
with the same spectral shape as model QG, but 10 times the in- 
tensity, and the dot-dashed contour is for the HMOl quasars only 
radiation field, scaled to the same hydrogen ionization rate as 
model QG. Each contour comprises two components, correspond- 
ing to net heating (low temperature component) and net cooling 
(high temperature component), respectively, which merge near 
the (density-dependent) thermal equilibrium temperature. 



5 THE RELATIVE IMPORTANCE OF 
DIFFERENT ELEMENTS 

For solar abundances, a large number of elements contribute 
to the radiative cooling rate. The black, solid curve in the 
left-hand panel of Figure |6] shows the normalized cooling 
rate (A/n^) as a function of temperature for a plasma 
in pure CIE and for solar abundances, re producing pre- 



vious w ork done on t hi s sub jec t fe.g.. | Cox fc Tucked 

1 19691: [Raymond et al." 'l976'; 'ShuU & van SteenberJ 

Il28y; iGaet z fc Salpeter 1983; Boehringer fc Henslerl 

1989;; S utherland fc Dopital Il993l; iLandi fc Landinil Il999l: 



Bcni amin et al.ll200ll ; ICnat fc Sternberell2007l : ISmith et all 
2008). The thin, coloured, solid curves show the contri- 
butions from individual elements, while the black, dashed 
curve shows the total cooling rate due to H and He. For 
T ^ 10^ K Bremsstrahlung (i.e., H & He) dominates, but at 
lower temperatures line radiation is most important. Going 
down from 10^ K to 10* K the cooling rate is successively 
dominated by iron, neon, oxygen, carbon, and hydrogen. 

If the cooling rate is dominated by a single element, as is 
for example the case for oxygen at T ~ 2 x 10^ K and for iron 
at T ~ 10^ K, then the total cooling rate will be sensitive 
to the relative abundances of those elements. For instance, 
since [O/Fe] is observed to vary with e nvironment by factor s 
of two or so at a fixed metallicity (e.g.. lShetrone et al.ll2003l ). 
we can expect similar variations in the cooling rates. 

We next turn on the HMOl ionizing background and 
illustrate the results for nu ~ lO"** cm~'^ and z = 3 in the 
right-hand panel of Figure [6] (recall that the left-hand panel 
is independent of both density and redshift). For this figure 
we have excluded Compton cooling off the CMB to isolate 
the impact of the ionizing radiation and because the former 
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Figure 6. Normalized cooling rates as a function of temperature for solar abundances, assuming either CIE {left-hand panel) or photo- 
ionization equilibrium for njj = 10~* cm"^ and an optically thin gas exposed to the z = 3 HMOl model for the meta-galactic UV/X-ray 
background from quasars and galaxies {right-hand panel). Note that normalized cooling rates are independent of the density for CIE, 
but not for photo-ionization equilibrium. The black, solid curve indicates the total cooling rate and the thin, coloured, solid curves show 
the contributions from individual elements. The black, dashed curve shows the contribution from H and He. 



is only important at high redshift. We can see that photo- 
ionization affects some elements more than others. As we 
have seen before, the effect is stronger for lower tempera- 
tures. Although we show the results for only a single den- 
sity here, we note that the importance of photo-ionization 
increases with decreasing density. 

Comparing the cooling rates including photo-ionization 
(right-hand panel of figure [6| to those for CIE (left-hand 
panel of figure O shows that photo-ionization increases the 
relative importance of oxygen and decreases that of carbon, 
helium, and especially hydrogen. It is also clear that many 
of the peaks of the various elements shift to lower tempera- 
tures when an ionizing radiation field is present. This shift 
occurs because a photo-ionized gas is overionized for its tem- 
perature compared to a coUisionally ionized plasma. If the 
ion fractions peak at lower temperatures, then so will the 
cooling rates due to coUisional excitation of those ions. 

This last figure illustrates the central result of this work: 
photo-ionization changes both the total cooling rates and 
the relative importance of individual elements. For a more 
complete visualization of this point, we kindly refer the in- 
terested reader to our web sitqij, where we host a number 
of videos, plots, and the tables themselves for download. 



6 DISCUSSION 

Radiative cooling is an essential ingredient of hydrodynam- 
ical models of a wide range of astrophysical objects, rang- 
ing from the IGM to (proto-)galaxies and molecular clouds. 
While numerical simulations of objects with a primordial 
composition often compute non-equilibrium radiative cool- 
ing rates explicitly and sometimes even include the effect 
of ionizing background radiation, the treatment of cooling 
of chemically enriched material is typically much more ap- 
proximate. For example, simulations of galaxy formation 



http : //www . strw . leidenuniv . nl/WSS08/ . 



typically either ignore metal-line cooling altogether or in- 
clude it assuming pure CIE. In addition, the abundances of 
all heavy elements are ty pically scaled by the same fa ctor 
(the metallicity) (but see iMartinez-Serrano et al.ll2008l and 
iMaio et al.l 120071 for recent exceptions) . In this simplified 
treatment metal-line cooling depends only on temperature 
and metallicity, allowing straightforward interpolation from 
pre-computed two-dimensional tables. 

We have used CLOUDY to investigate the effects of heavy 
elements and ionizing radiation on the radiative cooling of 
gas with properties characteristic of (proto-) galaxies and the 
IGM, i.e., optically thin gas with densities nn ^ 1 cm^'^ and 
temperatures r>10''K, assuming ionization equilibrium. 
We presented a method to incorporate radiative cooling on 
an element-by-element basis including photo-ionization by 
an evolving UV/X-ray background, using precomputed ta- 
bles, which for heavy elements are functions of density, tem- 
perature, and redshift and for H&He (which must be con- 
sidered together because they are important contributors to 
the free electron density) depend additionally on the He/H 
ratio. Using the 11 elements H, He, C, N, O, Ne, Mg, Si, S, 
Ca, and Fe, the redshift z = median absolute errors in the 
net cooling rate range from 0.33%, &t Z = O.IZq to 6.1% 
for the extreme metallicity Z = WZq, and the errors are 
smaller for higher redshifts. 

The tables as well as some scripts that illustrate how 
to use them are available from the following web site: 
http://www.strw.leidenuniv.nl/WSS08/. We also include 
tables for solar relative abundances which can be used if 
metallicity, but not the abundances of individual elements 
are known, as in equation ((5)1. This web site also contains a 
number of videos that may be helpful to gain intuition on 
the importance of various parameters on the cooling rates. 

We confirmed that, assuming CIE, heavy ele- 
ments greatly enhance the cooling rates for metallicities 
Z > 10~^ Zq and temperatures T < lO'^ K. We demonstrated 
that this remains true in the presence of photo-ionization by 
the meta-galactic UV/X-ray background. 

The background radiation removes electrons that would 
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otherwise be coUisionally excited, thus reducing the coohng 
rates. The effect is stronger for higher ionization parameters 
(i.e., higher radiation intensities or lower densities) and if the 
spectral shape of the radiation field is harder. Considering 
only the meta-galactic radiation field, which provides a lower 
limit to the intensity of the radiation to which optically thin 
gas may be exposed, the reduction of the metal-line cooling 
rates becomes important below 10® K for ionization param- 
eters U > 10~^ and below 10^ K for [/ > 10"^ (note that for 
the HMOl background U = 9 x W"'' /uh and 2 x 10"^/nH 
at z — and z = 3, respectively). 

As an example of the potential importance of includ- 
ing the effects of both photo-ionization and heavy elements, 
we considered the so-called warm-hot intergalactic medium 
(WHIM), which is thought to contain a large fraction of 
the baryons at redshifts z < 1. We demonstrated that the 
overdensities for which gas at typical WHIM metallicities 



(Z ~ 10 ^ Zq) and temperatures (T 



10' 



lO'^K) can 



cool within a Hubble time, can shift by an order of magni- 
tude depending on whether photo-ionization and metal-line 
cooling are taken into account. Hence, photo-ionization of 
heavy elements may have important consequences for pre- 
dictions of the amount of matter contained in this elusive 
gas phase. 

Because chemical enrichment happens in a number 
of stages, involving a number of processes with different 
timescales, the relative abundances of the heavy elements 
varies with redshift and environment by factors of a few. 
Hence, computing cooling rates on an element-by-element 
basis rather than scaling all elements by the metallicity, will 
change the cooling rates by factors of a few. The difference 
is therefore typically somewhat smaller than the effect of 
neglecting metals or photo-ionization altogether, but still 
highly significant. 

While it was known that different elements dominate 
the cooling for different temperatures in CIE, we showed 
that photo-ionization both shifts the peaks due to individ- 
ual elements to smaller temperatures and reduces their am- 
plitude. Note that since photo-ionization over-ionizes the 
gas, this effect is similar (but not equivalent to) that found 
in no n-equilibrium calculations without ionizing radiatio n 
(e.g., ISutherland fc Dopital [l993 : ICnat fc Sternberg||200'i1 ). 
Because the importance of photo-ionization depends on the 
ionization parameter, the relative contributions of individual 
elements exposed to a fixed ionizing radiation field depends 
also on the gas density. 

Would dropping our assumption of ionization equilib- 
rium have a large effect on the cooling rates? Ionizing ra- 
diation results in a plasma that is overionized relative to 
its temperature. Its effect is therefore similar to that of 
non-equilibrium ionization following rapid cooling (i.e., if 
the cooling time is shorter than the recom bination times 
of the ions dominating the cooling, see e.g.. iKafatod 119731 : 
IShapiro fc Moorelll976l ). We therefore anticipate that the ef- 
fect of non-equilibrium ionization will be much smaller for 
our cooling rates than for those that assume CIE. 

The assumptions that the gas is optically thin and ex- 
posed only to the meta-galactic background radiation are 
likely to be more important than the assumption of ion- 
ization equilibrium, particularly since non-equilibrium col- 
lisional cooling rates only differ from those assuming CIE 
by factors of a few or less (e.g., ISchmutzler fc Tscharnuteii 



1 19931 : [Sutherland fc Dopital [l993 : ICnat fc Sternberg! |2007| ). 
For column densities A^hi > 10^^ cm~^ self-shielding be- 
comes important and only part of the H-ionizing radiation 
will penetrate the gas cloud, which would particularly af- 
fect the cooling rates for T < 10^ K. At higher temperatures 
line cooling is dominated by heavier ions, which can only 
be ionized by higher energy photons and which therefore 
remain optically thin up to much higher column densities. 
This is because the photo-ionization cross sections of H and 
He drop rapidly with increasing frequency for energies ex- 
ceeding their ioinzation potentials. Moreover, for T ^ 10* K 
hydrogen is coUisionally ionized to a high degree and con- 
sequently the optical depth for ionizing radiation will be 
significantly reduced. 

It is, however, far from clear that high column den- 
sities would reduce the effect of radiation. For self-shielded 
clouds the cooling radiation may itself be trapped, providing 
a source of ioniz ing radiation even in the absence of an ex- 
ternal one (e.g., IShapiro fc Moorg|l976l : ICnat fc Sternberg 
[2007). Moreover, gas clouds with columns that exceed 
10^^ cm~^ are on average expected to be sufficiently close 
to a galaxy that local sources of ionizing radiation domina te 
over the background (|Schavell200^ : lMiralda-Escudell2005l ). 

Ultimately, these issues can only be resolved if non- 
equilibrium cooling rates are computed including radiative 
transfer and if the locations of all relevant sources of ionizing 
radiation are known. It will be some time before it is fea- 
sible to carry out such a calculation in, say, a cosmological 
hydrodynamical simulation. In the mean time, we believe 
that our element-by-element calculation of the equilibrium 
cooling rates for an optically thin gas exposed to the CMB 
and an evolving UV/X-ray background provides a marked 
improvement over earlier treatments. In future publications 
we will present cosmological, hydrodynamical simulations 
using these cooling rates. 
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